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■^J- , Abstract 

. If baryons couple only with w-mesons, one found the baryon spec- 

tral function may be negative. We show this unacceptable result is 
caused by the A^fc^-terms in the o;-meson propagator. Their contri- 
bution may not vanish in approximate calculations which violate the 
baryon current conservation. A rule is suggested, by which the calcu- 
lated baryon spectral function is well behaved. 
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Krein, Nielsen, Puff and Wilets (KNPW) [1] pointed out not long ago 
that in a self-consistent (SC) Hartree-Fock (HF) calculation of the renor- 
malized baryon propagator for the case of zero-density the spectral function 
Ar(k) can be negative for some real values of k, if baryons couple only with 
o;-mesons. They emphasized that this is unacceptable. The spectral repre- 
sentation they considered is of the form 

Since Ar(k) represents the probability that a state of mass \k\ is created, it 
must be non-negative. They suggested that it might be due to the inade- 
quacy of the HF approximation or the inconsistency of the theory. In their 
calculation they have neglected all the terms proportional to k^k v in the u- 
meson propagator on the basis of the baryon current conservation implied 
by the model for a rigorous calculation. Though this is a generally accepted 
approximation [2] and indeed, such terms need not be taken into account if 
the baryon current conserves, their contribution in the SCHF approximation 
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is not zero and has to be studied. It also indicates that the SCHF approxi- 
mation does not preserve the baryon current conservation. The aim of this 
letter is to show that the baryon spectral function becomes non-negative and 
is well behaved if an adequate part of the /c^/c^-terms is included. 

For a system of baryons coupling only with cj-mesons, the Lagrangian 
density has the form 

— 11 — 

L = -^(7,A + M)ip - -Fp,Fp, - -m 2 v A^ + ig^^A^ + L C tc, (2) 

where F^ u = d^A u - d v A^ = d/dx^, x^ = (x,ix ), x 2 = x^x^ = x 2 - 
Xq with x = t and CTC means the counterterm correction introduced for 
the purpose of renormalization. The baryon and u;-meson propagators are 
defined as 

G a p(x = Xl - x 2 ) = (T[V^i)V^2)]) = / jj^e^'Gapik), (3) 

D, u {x) = (T[A fl (xi)A v (x 2 )]) = J JiLe^'D^k), (4) 

where (0) = (^o|0|^o) an d l^o) denotes the exact ground state. The dressed 
HF scheme is illustrated in Fig. 1, where G (D) denotes an appropriate ex- 
pression chosen for the calculation of the baryon (u;-meson) propagator in the 
self-energy. As is wellknown, the relevant Dyson-Schwinger (DS) equations 
can be written as follows: 
(a) for baryon 



G{k) = G°{k) + G°(k)X(k)G(k) = -[7^ -iM + S(fc)]- 1 , (5) 
d T q 



£(*) = 9l J S 1 ^' 9, k - q)G(q)D vX (k - q) lx + Z CTC (k); (6) 



(b) for w-meson 

D, u {k) = D%{k) + D° m (k)U v x(k)Dxu(k), (7) 



%x(k) = -9l J T^I 7 ^* + q, k, q)G{k + q) lx G{q)i 



d T q 

where a caret indicates that the quantity is not yet renormalized, r = 4 — 
8 (8 — > + ) and in Eq. (5) the Feynman prescription M — * M — ie is 
understood. Since the baryon current is conserved, we have 

U v x(k) = (5 vX -^)U v (k), (9) 



2 



which gives U v (k) = |^n w (A;). Eq. (9) shows that the renormalized II^ A 
can be obtained from 

U v (k) = fl v {k) + UcTc(k). (10) 
The renormalized counterterms will be fixed by the following conditions 

dYi(k) 

S(A;)| 7M fc M=iMt = 0; #(7 k^) ^^= iM ^ = °' ( lla ) 



U v (k)\ k2=0 = 0; —^l\ k2=0 = 0, (lib) 



dU v (k) 
dk 2 

where M t is the true baryon mass. Following [3], we shall use the following 
notations to distinguish different approximations: scheme P, G = — [7„fc^ — 
iM t }-\ ~D Tj \ = D® x ; scheme BP, G — G, D rjX = L>° A and scheme FSC, 
G = G, D n x = D v \. It is seen that in scheme BP Eq. (5) has to be solved 
self-consistently [1, 4-6], while in scheme FSC [7, 3] one should consider a self- 
consistent coupled set of renormalized DS equations (5-10). As emphasized 
in [7], the study of an ansatz consistent with the Ward-Takahashi identity 
may provide a way to preserve the baryon current conservation. Though it 
does not hold in the SCHF approximation, we shall assume that Eq. (9) is 
valid approximately. We note that D^ u can be written as [8, 9] 

D^^-^^-i^^, (12) 

A v (k) = A° v (k) + A° v (k)U v (k)A v (k) = -i[k 2 +m 2 v + iU v (k) - ie]~\ (13) 

where 5m 2 is the mass counterterm for the cu-meson. As is wellknown [10], a 
properly normalized spectral representation for the baryon propagator equiv- 
alent to Eq. (1) can be written in the form: 



,-, + 7^ 2 l^a(-a 2 J Z 2 )+iM t p(-a 2 } Z 2 ) 
° {k) = -* VTW=* - J d ° k* + **-ie ' (14) 



where m\ = M t + m v is the threshold of the continuum spectrum, 

00 

Z 2 + J da 2 a(-a 2 ,Z 2 ) = 1 (15) 
and the spectral functions a and f3 should have the following properties: 
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(a) , they are all real, 

(b) .a(-<r 2 ,Z 2 )>0, (16) 

(c) . aa{-u\ Z 2 ) - M t (3(-a 2 , Z 2 ) > 0. 

For the cu-meson propagator we have 



AJk) = Z v Z! i Jdm 2 p ^~ m ' v \ (17a) 

" k 2 + m 2 — ie J k 2 + m 2 — ie 
th 



Z v + J dm 2 p v (-m 2 , Z v ) = 1. (17b) 

th 

In Eq. (17) m„ denotes the true cu-meson mass and th = (2M t ) 2 the threshold 
of the continuum. Since G(k) and G(k) as well as A v (k) and A v (k) are 
normalized differently, we should write G(k) = ZG(k) and A v (k) = Y v A v {k). 
It is easily seen that Z 2 = ZZ t and Z v = Y V R V , where (—Z t ) and (—iR v ) 
are the residues of G(k) and A„(&;) at the poles 7 M /c M = iM t and k 2 = — m 2 , 
respectively. In the case of zero-density one may write 

£(fc) = 7 M A; M a(A; 2 ) - iMb(k 2 ). (18) 

We note that under the on-shell condition (11a) one has Z t — 1 and M = M t 
[6]. Neglecting all the k fJi k u terms in Eq. (12) and substituting D^ u = 
S^ U A V = 5 l j LV Y~ 1 A v as well as G = Z~ X G in Eqs. (6) and (8), an explicit 
coupled set of renormalized equations for the determination of (a, b), (a,/3), 
U v and p v has been given in [3]. Since the expressions for the rest of the 
equations considered here are the same as in [3], for space saving and for 
convenience of discussion only a part of them will be written down as follows: 

a(k 2 ) = a s (k 2 ) + a A (k 2 ); b(k 2 ) = b s (k 2 ) + b A (k 2 ) (19) 



a s(k 2 ) = ^jj da 2 dm 2 J dxf a (-a 2 )h v (-m 2 )(l-x) In ^ t } + c s , 
oo ^ ' 

(19a) 

bs(k 2 ) = {^Jj da 2 dm 2 J dxf p (-a 2 )h v (-m 2 ) In £ « ] + c 5 , (19b) 
oo ^ ' 

/ 7 (-a 2 ) = 5{a 2 - M 2 ) + 9(a 2 - m 2 ) 7 (-a 2 ), ( 7 = a or /?) (19c) 
h v (-m 2 ) = R v 5(m 2 - m 2 v ) + 6(m 2 - th)p v (-m 2 ), (19d) 
K 2 (k 2 ) = K 2 (x, k 2 , a 2 , m 2 ) = x(l - x)k 2 + xa 2 + (1 - x)m 2 , (19e) 
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where 7(-o" 2 ) = Z 1 ^(-a 2 ,Z 2 ) and p v (-m 2 ) = Y v 1 p v (-m 2 , Z v ) . In addi- 
tion, we have 



/, 2 x 1. l + a(k 2 ) „ /l2 . 1 1 + 6(A; 2 ) . . 

a(k 2 ) = -/m-^i; ^(A; 2 ) = -/m-^i, (20a) 

L>(/c 2 ) = [1 + a(A; 2 )] 2 A: 2 + [1 + 6(fc 2 )] 2 Mf. (20b) 

In Eq. (19) the A-terms are contributed by the second A^A^-term in the 
round brackets of Eq. (12). They are given by 

P P OO P 1 /' 1_ — t' 

a A (k 2 ) = -r&p JJ ^ da 2 dm 2 J ^dx J ^ dyf a (-a 2 )h v (-m 2 )x 
{(1 + 3x) In ^1 - * 2 (1 - x)[^ + + c A 



b{±(k 2 ) = -^JJ da 2 dm 2 J dx J dyfp(-a 2 )h v (-m 
ro ln ^ 2 (-V 2 ) I r 2 r k 2 , m| n 



(21b) 

L 2 (fc 2 ) ^ X ^L 2 (fc 2 ) ^ L 2 (-M 2 )J ^ 

L 2 (A; 2 ) = L 2 (x, y, A; 2 , a 2 , m 2 ) = x(l - x)k 2 + :ra 2 + ?/m 2 , (21c) 

where in Eqs. (19) and (21) we have c K = 2M 2 [a' K (-M 2 ) - b' K (-M 2 )}, k = 5 
or A, and f'{k 2 ) = df(k 2 )/dk 2 . 

As shown in Eq. (20), in order to find a and (3 we need to know 
the imaginary parts of a and b. Eqs. (19) and (21) show that a and 
b will be complex if either K 2 (k 2 ) or L 2 (k 2 ) becomes negative. Define 
Q K = 9 K (x,a 2 ,m 2 ,k 2 ) = if K 2 (k 2 ) > and 6 K = 1 if K 2 (k 2 ) < 0. We 
may write lnK 2 (A; 2 ) = In \K 2 {k 2 )\ ± i9 K (n + 2rar). Even if we set n — 0, 
ln.fr 2 is still indeterminate in case K — 1. Let the sign before iir be denoted 
by slnfT 2 . Obviously the same remark also applies to In L 2 (k 2 ). We shall 
make the stipulation that sln.fr 2 and slnL 2 are to be so chosen that Eq. 
(16) should be satisfied. Since a and b are derived from the same righthand 
side of Eq. (6), their slni^ 2 (slnL 2 ) should be the same. According to [3], 
the above coupled set of equations can be solved easily by the method of 
iteration. We have solved it for schemes BP and FSC. In the following we 
shall choose M t = 4.7585/m^ 1 , m v /M t = 0.83387 and g 2 v = g 2 /8n 2 = 1.3685. 
If we neglect all the A^Av-term in Eq. (12), it is impossible to make a(A; 2 ) 
non- negative by an appropriate choice of slnfT 2 and slnL 2 . This is unac- 
ceptable and confirms the result of KNPW. As is wellknown, the last term in 
Eq. (12) is not renormalizable and the model becomes renormalizable only, 
because the tu-meson couples to the conserved baryon current so that this 
term may be neglected. Thus we have calculated the additional contribution 
of the second term in the round brackets of Eq. (12) according to Eq. (21). 
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The numerical results for sin K 2 = — and sin L 2 = + are shown in Fig. 2a. 
They are the best results we can find and satisfy all the requirements of Eq. 
(16). However, they are still physically unacceptable, because the threshold 
is wrong. It should be —(M t + m v ) 2 rather than —M 2 . This implies that the 
results contain a contribution from an additional zero-mass meson, which 
does not exist in our model. Clearly the latter is caused by the fact that 
A v contains an additional zero-mass pole, because k 2 = k 2 — k$ = does 
not ensure that k^ and k v are both zero. Since Eq. (12) is formally rigorous 
and D^ u should contain no zero-mass pole, the renormalized contribution of 
the latter must be cancelled by a finite part arising from the last term. In- 
deed, this can also be seen from the zero-order approximation to D^, which 
reads 



DIM = - -ir)AS(*) - n£4 = (<W + -Ar)AS(*), (22) 



where A°(fc) = — i[k 2 + ml — ie] -1 . Thus, if we neglect the last term in Eq. 
(12), D Q will contain a false zero-mass pole. Eq. (22) indicates that this 
last term makes the pole structure of D Q correct, though unrenormalizable. 
Our task is then to eliminate the contribution of the zero-mass meson from 
Eq. (21). This can be done unambiguously. We note that the integrals over 
the parameter y in Eq. (21) can be performed exactly and in an elementary 
way. The expressions are, however, too lengthy to write down here. Let the 
curved brackets in Eqs. (21a and b) be denoted by Y c (y) with c = a and b, 
respectively. We may express the integration over y as 



Set L 2 (y) = L 2 (x,y,k 2 ,a 2 ,m 2 ). Eq. (19e) and (21c) show 

L 2 (y = 0) = K 2 (x,k 2 ,a 2 ,0); L 2 (y = 1 - x) = K 2 (x, k 2 , a 2 , m 2 ). (24) 

Form Eqs. (21), (23) and (24) one easily finds that the contributions to the 
wrong threshold come exclusively from Z c (0). Substituting Eq. (23) in Eq. 
(21), discarding the contribution of Z c (0) and using Eq. (19), we have again 
solved the coupled set of equations. Fig. 2b depicts our numerical results for 
a(k 2 ) and f3(k 2 ) with the choice sln.fr 2 = — and slnL 2 = — . They satisfy 
all the requirements of Eq. (16) and have a correct threshold. Thus they 
are now acceptable. As displayed in Fig. 2, the effect of self-consistency on 
the baryon spectral function is perceptible, though the results of schemes 
BP and FSC are almost the same. The latter asserts, in agreement with [3], 
that there is no need to include the meson propagator in the self-consistency 



k,,k, 



'V 



'V 




(23) 
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requirement. Clearly the fulfillment of Eq. (16) also implies that A r (k) will 
be non-negative. Our results show if in an approximation the contribution of 
the renormalizable k^k^-term in D^ u is not zero, it must be taken account of, 
otherwise the result may be qualitatively wrong. According to our above re- 
sults the rule for the calculation of this additional contribution may be stated 
as follows: for c = a or 6, Z c (0) and the last term in Eq. (12) should both be 
discarded, because the contributions of the former and of a finite part from 
the latter cancel each other, while for the sake of renormalization the infinite 
nonrenormalizable part need not be considered on the basis of baryon current 
conservation. The latter can be understood as follows. Let A = A + A r + 
be a term calculated in the SCHF approximation, where the subscripts 0, r 
and oo denote the contributions from the first, second (renormalizable) and 
third (nonrenormalizable) term in Eq. (12), respectively. Since the model 
possesses the baryon current conservation, if A violates it, there must exist 
another term B in the theory such that A + B preserves it. Thus, A r + B r 
and Aqo + Boo can be neglected. Since the infinity is an indefinite quantity 
and it always contains an undefined finite part (indicated by /), we may 
write Aao + B^ = A'^ + B'^ + Af + Bf. Thus, we are allowed to throw away 
A'qq + B'qq, if an adequate rule is prescribed to retain Af + Bf. If only A is 
calculated and it can already yield a quite good approximation, clearly A r 
and Af must be taken into account except when they are very small. 

Redmond [11] suggested that one may use the spectral representation to 
eliminate the ghost poles. Here we would like to point out that there is a 
simple and direct way to derive a ghost-free representation for the Dyson 
solution by the spectral representation. From Eqs. (5), (14) and 

G(k) = -D- 1 ^^ (1 + a) + iM t (1 + ft)] = Z~ x G(k), (25) 

where D is given by Eq. (20b), we easily find 

l + a = F{k 2 )- l F a {k 2 )] 1 + 6 = F(k 2 y l Fp(k 2 ), (26a) 

F(k 2 ) = k 2 F a (k 2 ) 2 + M 2 F^(k 2 ) 2 (26b) 

with F 1 (k 2 ) = J °° da 2 f^(-a 2 )(k 2 + a 2 - ie)~ l (7 = a or /?). It is seen that if 
a(k 2 ) and f3(k 2 ) are known, so are F a {k 2 ) and Fp(k 2 ). Let us use a subscript 
R to indicate the results obtained from Eq. (26). Since a R and 6^ are direct 
consequences of the spectral representation G{k), G R (k) deduced from these 
<2r and 6# by means of Eq. (5) satisfies Eq.(25) rigorously and is thus 
ghost-free, because G(k) is well behaved and possesses the Herglotz property 
[12]. Further one observes that G R (k) obviously satisfies the following Dyson 
equation: 

G R (k) = G°(k) + G°(k)T, R (k)G R (k), (27) 
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with G°{k) = -{-i^k^-iMt}- 1 and T, R {k) = 7 ^k^a R (k 2 ) - iM t b R (k 2 ). Let us 
insert a and (3 found from Eq.(20) with a = 0,5 + a a and b = 65 + 6a calculated 
according to Eqs.(19) and (21) in F a and Fp. From Eq. (26) we obtain a R 
and b R . We note that the set (a R , b R ) is different from (a, b), though the 
former is deduced from the latter. In fact, if we substitute (a, b) into Eq. 
(25), the latter no longer holds, because G(k) found in this way possesses 
ghost poles, which are —100.4468 ± i55.2553/m -2 . The reason is clear. To 
calculate a and b we need renormalization, which, however, does not ensure 
that Eq. (25) holds. In Fig. 3 we have compared {a R ,b R ) with (a,b). They 
have different asymptotic behavior and differ widely, though the same a and 
P can be derived from them by Eq. (20). Our numerical result for p v (k 2 ) is 
well behaved and is qualitatively similar to p v (k 2 ) found for the a — u model 
in [3]. Using Eqs. (13, 17) and the relation A v (k 2 ) = Y~ 1 A v (k 2 ), we can 
also derive an analytic expression for U v (k 2 ), which will yield a ghost-free 
A„ (k 2 ) by means of Eq. (13). However, for lack of spacing both p v and 
11^ will not be discussed here. As shown in Refs. [1, 3, 4, 7], if in addition 
to the u;-meson, other mesons, for instance n, a and chiral n — a (linear 
cr-model), are considered, the baryon spectral functions in the self-consistent 
HF approximation are regular, even though the /c M /c„-terms in Eq. (12) are 
neglected. The question what their effects are in these more complicated 
cases is being studied. 

The work is supported in part by the National Natural Science Foundation 
of China and the Foundation of Chinese Education Ministry. 
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Fig. 1: Diagrammatic representation of the dressed HF approximation, 
(a) the baryon propagator; (b) the oo-meson propagator. 
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Fig. 2: The baryon spectral functions a(k ) and (3(k ). 
(a) with wrong threshold; (b) the correct final result. 
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Fig. 3: Comparison betwwen (a(k 2 ), b(k 2 ) and (a R (k 2 ), b R (k 2 )). 

The subscript V for the real part, while the T for the imaginary part. 



